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A measurement is presented of the inelastic proton-proton cross section at a centre- 
of-mass energy of ^/s = 7TeV. Using the CMS detector at the LHC, the inelastic 
cross section is measured through two independent methods based on information 
from (i) forward calorimetry (for pseudorapidity 3 < \r]\ < 5), in collisions where at 
least one proton loses more than 5 x 10^^ of its longitudinal momentum, and (ii) the 
central tracker ( < 2.4), in collisions containing an interaction vertex with more 
than 1, 2, or 3 tracks with transverse momenta pj > 200MeV/c. The measurements 
cover a large fraction of the inelastic cross section for particle production over about 
9 units of pseudorapidity and down to small transverse momenta. The results are 
compared with those of other experiments, and with models used to describe high- 
energy hadronic interactions. 
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1 Introduction 

Total hadronic cross sections, as well as their major subdivisioris ir\to elastic, inelastic diffractive 
and inelastic non-diffractive contributions, comprise fundamental quantities that have been 
studied in high-energy particle, nuclear, and cosmic-rays physics over the past 60 years, in 
experiments covering many orders of magnitude in centre-of-mass energy ||l]-(5l. 

The bulk of the total cross section in proton-proton (pp) hadronic interactions cannot be cal- 
culated through perturbative quantum chromodynamics, but phenomenological approaches 
based on fundamental principles of quantum mechanics, such as unitarity and analyticity, 
can be used to accommodate the experimental results (e.g. Ref. (SJ, and references therein). 
Although phenomenological models of cross sections at low centre-of-mass energies {\/s < 
100 GeV) provide a rather precise description of the data, there are large uncertainties in ex- 
trapolating to the energy range of the Large Hadron Collider (LHC). The measured inelastic pp 
cross section ((Ji^ei) serves as an input to these phenomenological models, and provides basic 
information needed for tuning hadronic Monte Carlo (MC) generators. The values of t7inei are 
also used to estimate the number of pp interactions as a function of luminosity at colliders, 
and are relevant to studies of high-energy cosmic rays [Zl and to the characterization of global 
properties of heavy-ion collisions, especially in the context of the Glauber model [SI. 

This Letter presents a measurement of the inelastic pp cross section at = 7 TeV, using data 
collected with the Compact Muon Solenoid (CMS) detector at the LHC. The analysis is based 
mostly on the central silicon tracker and the forward hadron calorimeters (HF) of the CMS 
apparatus. The combination of these two detectors provides sensitivity to a large part of the in- 
elastic cross section, including central diffractive production, where particles can be produced 
at small values of pseudorapidity. 

The measurement using the HF calorimeters covers a region of phase space corresponding to 
values of fractional momentum loss of the scattered proton of ^ = {MxC^)^/s > 5 x 10^^, 
equivalent to Mx > 16GeV/c^, where Mx is defined as the larger mass of the two dissociated 
proton systems in the final state. This coverage is the same as that used in recent publications 
by the ATLAS f3] and the ALICE |5| Collaborations. 

2 Experimental apparatus 

A detailed description of the CMS apparatus can be found in Ref. 0, and the features most rel- 
evant to the present analysis are sketched below. The CMS detector comprises a 6 m diameter, 
13 m long, 3.8 T solenoid magnet, with a combined silicon pixel and strip tracker covering the 
region \rj\ < 2.5, a lead-tungstate electromagnetic calorimeter and a brass/ scintillator hadronic 
calorimeter covering the region \ f]\ < 3.0; these detectors are contained within the volume of 
the magnetic field. The pseudorapidity is defined as = — In [tan [9/2)], where 9 is the po- 
lar angle of any particle with respect to the anticlockwise circulating beam. Several layers of 
muon chambers (drift-tube, resistive-plate and cathode-strip chambers) form the outer part of 
the detector. The charged-particle resolution of the central tracker for a transverse momentum 
of Pt = 1 GeV/c is between 0.7% at |?7 1 = and 2% at |f/ 1 = 2.5 |19|. 

On each side of the detector, at 3.0 < \ri\ < 5.2, reside the hadron forward calorimeters (HF), 
each composed of 18 iron wedges, with embedded quartz fibres running along the beam direc- 
tion. Each wedge is subdivided into 13 f/-segments, called towers. 

The beam-sensitive "pick-up" detectors, consisting of two pairs of button electrodes located at 
±175m from the centre of the detector, provide almost 100% detection efficiency and accurate 
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3 Estimating the inelastic cross section using the IHF calorimeters 



timing of proton bunches at CMS. The luminosity is calculated with a precision of ±4% by 
relating the energy depositions in the HF calorimeters to the results obtained during dedicated 
van der Meer scans 1121 HH- 

3 Estimating the inelastic cross section using the HF calorimeters 

In this method, the inelastic pp cross section is measured by counting the number of events that 
deposit at least 5 GeV of energy in either of the two HF calorimeters. The threshold Ehf > 5 GeV 
is set to minimize the effect of detector noise on the efficiency of selecting pp collisions. 

3.1 Event selection and analysis 

The analysis is performed using data collected in low-luminosity runs with an average of 0.007 
to 0.11 collisions per bunch crossing. The events are collected using three triggers: (i) a coin- 
cidence trigger that requires the presence of two colliding bunches, used to select an unbiased 
sample of pp events, (ii) a single-bunch trigger, requiring the presence of just one unpaired 
bunch, used to estimate beam-induced backgrounds, and (iii) a random "empty" trigger, re- 
quiring absence of both beams, which is used to estimate detector noise. All these triggers are 
formed from information provided by the beam pick-up detectors. 

The analysis is based on counting the number of pp collisions with £hf > 5 GeV in either of 
the two HF calorimeters. The cross section is evaluated in terms of the variable ^, which is 
defined through MC studies as follows. For each MC event, generator-level information is 
used to order final-state particles in r] and to find the largest gap in r] between two consecutive 
particles. This "central" gap is used to separate all particles into two groups, by assigning each 
particle, according to its rj position relative to that gap, to system A or system B. Finally, the 
masses of system A and B are calculated, and the larger of the two is called Mx, while the 
smaller one My, thereby defining ^ = (Mxc^)^/s. In single-diffractive events, ^ corresponds 
to the fraction of momentum lost by the proton in a collision. 

The distributions in ^ values for £hf > 4 and > 5 GeV are shown in Fig.[l]for three Monte Carlo 
models: PYTHIA 6 (version 6.422) [12J, PYTHIA 8 (version 8.135) |13|, and Phojet (version 1.12- 
35) imillSl. As the plots illustrate, to maintain large detection efficiency, and to mitigate model- 
dependence, it appears adequate to restrict the range of ^ to values greater than 5 x 10^^. The 
measured values of CTj^ei are corrected using two quantities obtained through MC simulation: 
the selection efficiency e^, which is the fraction of pp interactions with ^ > 5 x 10^^ that is 
selected by requiring Ehf > 5 GeV, and the contamination which is the fraction of events 
that have Ehf > 5 GeV, but originate from ^ < 5 x 10^^. Table [l] gives the values of eg and 
estimated in the three Monte Carlo models. As the table shows, the criterion Ehf > 5 GeV 
selects a large fraction of events with ^>5xl0~^, with only a small contamination from events 
with ^ < 5 X 10^^ that characterize contributions originating from low-mass single-proton or 
double-proton fragmentation. 

3.2 Measurement of the inelastic cross section 

The analysis is performed using 9.2 million events, corresponding to an integrated lumi- 
nosity of 2.78 f/b^^, collected under the two-bunch coincidence condition, of which 2.1% have 
Ehf > 5 GeV. The fractions of Ehf > 5 GeV events selected by the single-bunch and empty 
triggers are, respectively, 0.30% and 0.32%, suggesting that most of the single-bunch events are 
from detector noise rather than beam-gas collisions. This is confirmed by the observation that. 
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Figure 1: The normalized ^ distributions for £hf > 4 and £hf > 5 GeV from MC simulation of 
inelastic pp collisions using (a) PYTHIA 6, (c) PYTHIA 8, and (e) Phojet, are shown for the full 
range of ^. The corresponding efficiencies are shown in (b), (d), and (f), respectively. The cutoff 
value of ^ used in this analysis of 5 x 10^^ is shown on the plots as a dashed vertical line. 
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3 Estimating the inelastic cross section using the IHF calorimeters 



Table 1: Values of efficiency {e^) and contamination (b^) for events with ^ > 5 x 10^^ using 
the selection criterion of Ehf > 5 GeV, obtained for three Monte Carlo models of hadronic 
production. 



Generator 


6^ (%) 


h (%) 


PYTHIA 6 


97.5 ± 0.6 


2.0 


PYTHIA 8 


99.3 ± 0.2 


2.0 


Phojet 


99.1 ± 0.2 


1.2 



in the single-bunch triggered sample, the number of events with at least one track is very small. 
For this reason, beam-gas contributions are considered negligible. 

The number of detected inelastic collisions (Nj^el) contained in the total number of coincident 
inelastic events (Ncoinc) is obtained as follows: 

Ninel — ^comc[(-^ coinc F empty) F empty (-^ :oinc F empty)]/ (1) 

where J-'empty arid J-'coinc correspond to the fractions of empty and coincidence triggers with 
Ehf > 5 GeV. The term Ncoinc-^empty(-^coinc — -^empty) represents the number of true collisions 
in Ncoinc-^^empty events. 

The value of Nmei has to be corrected for event pileup, i.e. the possibility that more than one 
collision with Ehf > 5 GeV occurs in the same trigger, but all such collisions are counted as just 
a single event. The number of collisions per trigger is assumed to follow Poisson statistics, for 
which the probability of / simultaneous collisions (/ = 1, 2, 3, . . .) is given by 

n^.A) = -^, (2) 

where A is the mean number of interactions with Ehf > 5 GeV, which depends on the instan- 
taneous luminosity (L). The fraction /pu of overlapping collisions, each with Ehf > 5 GeV, is 
computed as 

f __ ^ZiPji'^) _ 1-(1+A)e-^ 
/pu - E~iP(i,A) ~ l-e-'\ 

-J-i + 0{A^), (3) 

where A is evaluated from the fraction of detected interactions rint = Ninei/ Ncoinc ^ 

rini = E~iP(/,A) = 1-P(0,A) = l-e-\ 
A = -ln(l-rint). (4) 

The denominator in Eq. ^ assumes independent probabilities for detecting each of the simul- 
taneous collisions, which is a good approximation for Ehf > 5 GeV. 

Table |2] lists the values of A and /pu and their statistical uncertainties for different data runs. 
The accuracy on the correction factor /pu is limited mostly by the number of events in each run. 

The relationship used to evaluate the cross section for ^ > 5 x 10^^, taking account of correc- 
tions for pileup, efficiency, and contamination corresponds to: 

rr (^^^^^n-^^ Mnel(l - &g)(l +/pu) 

Ci^eK^ > 5 X 10 ) = ^^j^^^ . (5) 

where J Ldt is the integrated luminosity of the data sample. 
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Table 2: Mean number of collisions with Ehf > 5 GeV per coincidence trigger (A) and fraction 
of overlapping collisions (/pu) for the runs used in this analysis. 



Run No. A /] 



132601 (0.64 ± 0.01) % 0.0032 ± 0.0001 

132599 (0.78 ± 0.01) % 0.0039 ± 0.0001 

133877 (1.74 ± 0.02) % 0.0087 ± 0.0001 

133874 (3.34 ± 0.05) % 0.0166 ± 0.0002 

137027 (4.59 ± 0.17) % 0.0228 ± 0.0009 

135575 (8.41 ± 0.04) % 0.0415 ± 0.0002 

135175 (9.98 ± 0.05) % 0.0491 ± 0.0003 



3.3 Results and systematic uncertainties 

The value of (Tjnei for ^ > 5 x 10^^ is calculated by averaging the results obtained from Eq. ^ 
for the different pileup conditions of Table |2j The largest systematic uncertainty, besides the 
4% uncertainty on luminosity, is from the observed run-to-run changes in the parameters of 
Eq. (|5|, which we take as an additional systematic uncertainty of the measurement. The model 
dependence of the efficiency contributes ±1%, while the correction for the contamination 
from events below the ^ threshold is uncertain by ±0.5%. The exclusion of noisy HE towers in 
the calculation of HE energy changes the results by ±0.4%, a value that is taken as a system- 
atic uncertainty. Einally, lowering the value of the calorimeter threshold £hf from 5 to 4 GeV 
introduces a change of 0.2% in the final result. 

Table 3: List of systematic sources and their effects on the value of the inelastic cross section 
measured using HE calorimeters. The luminosity contributes an additional uncertainty of 4% 
to this measurement. 



Systematic source 


Uncertainty on cTinei 


Change in cT-nei 


Run-to-run variation 


±0.8 mb 


±1.3% 


Selection efficiency 


±0.6 mb 


±1.0% 


Contamination from ^<5xl0^^ 


±0.3 mb 


±0.5% 


HE tower exclusion 


±0.3 mb 


±0.4% 


HE energy threshold 


±0.1 mb 


±0.2% 


Total (in quadrature) 


±1.1 mb 


±1.8% 



Table|3]lists the individual systematic uncertainties, and their total impact, calculated by adding 
the separate contributions in quadrature. The inelastic pp cross section for events with ^ > 
5 X 10"^ is found to be: 

(7i^ei(^ > 5 X 10-6) = [60.2 ± 0.2 (Stat.) ± 1.1 (syst.) 

± 2.4 (lum.)] mb. 

This result is in agreement with equivalent measurements from the ATLAS Collaboration 
(7-ATLAS(^ > 5 X 10-6) = [60.3 ± 0.05 (stat.) ± 0.5 (syst.) ± 2.1 (lum.)] mb [|3J, and from the ALICE 
Collaboration chalice > 5 ^ IQ-^) = [62.12509 (syst.) ± 2.2 (lum.)] mb ||5l. The uncertainties 
on luminosity of the three measurements are highly correlated. 
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4 Estimating the inelastic cross section by counting event ver- 
tices 

A vertex-counting method is also used to measure the inelastic pp cross section. The method 
relies on the accuracy of the CMS tracking system and not upon any specific Monte Carlo 
simulation. This method assumes that the number (n) of inelastic pp interactions in a given 
bunch crossing follows the Poisson probability distribution of Eq. (|2|, where A is calculated 
from the product of the instantaneous luminosity for a bunch crossing and the total inelastic pp 
cross section: A = L ■ (Jjng\. The probability of having n inelastic pp interactions, each producing 
a vertex with >1, >2, or >3 charged particles with pj > 200MeV/c within \r]\ = 2.4, for n 
between and 8, is measured at different luminosities to evaluate t7inei from a fit of Eq. Q to 
the data. 

4.1 Event selection and method of analysis 

Inclusive samples of ~3 x 10^ two-electron candidate events, and ^1.5 x 10^ single-muon can- 
didate events, are selected for this analysis. The specific trigger requirements are not important, 
as long as their efficiencies do not depend on the number of pileup interactions. The "trigger- 
ing interaction", i.e. the process associated with the production of either the two electrons or 
the single muon, is not included in the vertex count, but is used just to sample unbiased pileup 
interactions, given by the additional vertices in the same bunch crossing. 

For each of the two data samples, the distributions in the fraction of events with to 8 pileup in- 
teractions are measured as a function of luminosity. A bin-by-bin correction is applied to these 
measurements to account for possible migrations and inefficiencies. The final distributions are 
then fitted to Eq. (|2|, to extract a common value of (Tjj^eh 

The bin-by-bin corrections, evaluated from full Monte Carlo simulation (PYTHIA 6) and recon- 
struction of events in the CMS detector, do not depend on any specific production model, but 
only on an accurate simulation of the CMS tracking system. The distributions of charged tracks 
in transverse momentum and in track multiplicity in MC events are reweighted to provide 
agreement with the data, as these two quantities influence the vertex reconstruction efficiency. 
Vertices with a large number of tracks are always well reconstructed, while vertices with less 
than 10 tracks suffer some degradation in reconstruction efficiency. Cross sections are measured 
for inclusive pp interactions with >1, >2, and >3 charged particles, with pj > 200MeV/c and 
< 2.4, where "charged particles" are referred to those with decay lengths ct > 1 cm. 

4.2 Vertex definition and reconstruction 

To be counted, a pileup interaction has to have a sufficient number of tracks to provide a vertex 
of good quality [16|. The vertex quality depends upon the number and characteristics of the 
individual tracks attributed to each vertex. 

There are two main reasons that lead to incorrect vertex reconstruction: (i) overlap with another 
vertex, i.e. the reconstruction program merges two vertices, and (ii) an insufficient number 
of tracks, or tracks too poorly measured to pass vertex-quality requirements. The vertexing 
algorithm is very efficient in distinguishing vertices that are further apart than 0.06 cm along the 
beam direction, and this analysis requires a minimum distance of 0.1 cm. Minimum distances 
of 0.06 cm and of 0.2 cm are used to check for any systematic effects from this requirement. The 
fraction of vertices lost from merging depends on luminosity is as large as 2% in these data, 
an effect that is well reproduced by MC simulation, and is therefore corrected. The second 
source of inefficiency in vertex reconstruction depends on the number of tracks per vertex. 



4.3 Results and systematic uncertainties 
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The reconstruction efficiency is about 40% for a 2-track vertex, and reaches 100% for 10-track 
vertices. This is accounted for in the calculation of efficiency. 

There are also two main reasons for generating secondary vertices that are not related specif- 
ically to cTinei: additional vertices generated through decays of long-lived particles, and false 
secondary vertices generated by splitting a single vertex into two distinct vertices. Misidenti- 
fied secondary vertices can often be rejected, as they have a much lower track multiplicity, and 
they are not necessarily positioned along the beam line. 

The correction of number of candidate vertices to the true number of pileup interactions is con- 
sidered as a function of luminosity. In particular, a 2-vertex event recorded at low luminosity 
is most likely to correspond to a true 2-vertex event, while a 2-vertex event recorded at high 
luminosity is most likely a 3 or 4-vertex event, in which 1 or 2 vertices are merged. We divide 
the data into 13 equal intervals of instantaneous bunch-crossing luminosity, from 0.05 x 10^^ 
to 0.7 X 10'^'^ cm^'^ s^^. To obtain the true pileup distribution in each luminosity interval, we 
proceed as follows: 

(i) Using Eq. the expected distribution of pileup interactions is calculated for the specific 
luminosity interval, assuming some trial value crl^^i for the inelastic cross section. 

(ii) The Monte Carlo simulation is reweighted to generate a pileup distribution matching the 
one calculated in step (i). Steps (i) and (ii) are repeated several times for different c^^^gf , until 
good agreement is reached between data and the reconstructed pileup distributions for MC 
events. 

(iii) The generated pileup distributions for inclusive interactions with >1, >2, and >3 tracks, 
each with pj > 200 MeV/c and \ < 2.4, is obtained from the reweighted Monte Carlo. 

(iv) The bin-by-bin corrections are computed using the ratio of reconstructed to generated 
Monte Carlo pileup distributions for >1, >2, and >3 tracks, yielding thereby the correction 
factors for each of these three inclusive sets of events. 

The corrected fractional distributions of events, for interactions with more than 1 track in data 
or in the MC, are compared in Fig. |2]as a function of the number of vertices (n) for the eleven 
lowest bins in instantaneous luminosity. 

4.3 Results and systematic uncertainties 

Figure |3] displays the data points from Fig.|2]as a function of the instantaneous luminosity, for 
events with n = to 8 pileup vertices. For each n, the values of the Poisson distribution given 
by Eq. (|2| for i = n are fitted as a function of A = L ■ (Jmei to the data, providing 9 estimates 
of the inelastic cross section. Their weighted average provides the final result shown in Fig. |4] 
(a). The error bars and the values of goodness of fit per degree of freedom (x^/NDOF) for each 
result are obtained from the individual Poisson fits of Fig.|3] Figure |4](b) shows the normalized 
values for each of the fits to Eq. 

The main source of systematic uncertainty is the 4% uncertainty on CMS luminosity, which 
leads to an uncertainty Aai^j^ = ±2.4 mb. The largest contribution arising from the method 
of analysis, AcTVtx = ±1.4 mb, is the uncertainty on the vertex-reconstruction efficiency, which 
is evaluated using a Monte Carlo simulation and a method based on data. This second tech- 
nique utilizes measured quantities such as the distribution of the longitudinal z position of 
the vertex and the distribution of the minimum distance between two vertices to evaluate the 
vertex-reconstruction efficiency. Other uncertainties linked to vertex selection are estimated 
by: (i) reducing the range used for accepting longitudinal positions of vertices from 20 to 10 cm 
relative to the centre of the CMS detector, (ii) modifying the vertex-quality requirements, (iii) 
changing the minimum distance between two vertices from Az < 0.1 cm to 0.06 cm and 0.2 cm. 
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Figure 2: Fraction of reconstructed events with more than 1 track, corrected for efficiency, mea- 
sured as a function of the number of vertices, in data (dots) and in Monte Carlo (histogram), 
for instantaneous bunch-crossing luminosities between 0.05 x 10^^ and 0.6 x 10'^'^ cm^'^ s^^. 

and (iv) changing the maximum transverse coordinate of the vertex from ±0.06 cm to ±0.05 
cm and ±0.08 cm. Several other possible sources of uncertainty have also been checked by: 
(i) performing the analysis on two different sets of data, (ii) changing the luminosity interval 
used in the fit by ±0.05 x lO^'^cm^-^s^^, and (iii) repeating the analysis without reweighting 
the track-multiplicity distributions in the MC. The uncertainty attributed to each systematic 
source is defined by the largest change in tTinei- The full list of the systematic sources is shown 
in Table |4j Adding all the uncertainties in quadrature yields a total systematic uncertainty on 
the method of Ac^syst = ±2.0 mb. 

The measured values of tTi^ei for inclusive interactions with >1, >2, and >3 charged particles 
with \r]\ < 2.4 and px > 200MeV/c, as well as their individual uncertainties, are listed in 
Table il 

5 Results and comparison with Monte Carlo models 

The two techniques presented to measure the inelastic pp cross section complement each other. 
The calorimeter-based method is very sensitive to events that produce forward energy de- 
position, and, in particular, small Mx values that comprise particle systems highly boosted 
along the beam line. However, the method is less sensitive to central diffractive dissocia- 
tion events, with particle production concentrated at small pseudorapidities. Conversely, the 
vertex-counting method is geared toward measurement of centrally-produced events, and is 
not optimal for events with particles produced mostly at large rj. The concurrent use of these 
two methods provides therefore almost complete coverage of all types of pp inelastic events, 
with particle production in the range of |?/ 1 < 5. 
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Figure 3: Fraction of pp events with n pileup vertices, for n = to 8, containing more than one 
charged particle, as a function of instantaneous bunch-crossing luminosity. The dashed lines 
are the fits described in the text. The data points are plotted at the mean of the differential 
distribution in each bin. 

Figure |5] compares the CMS results with the measurements presented by the TOTEM HI, the 
ATLAS [3J and the ALICE [5J collaborations, as well as with predictions of two groups of Monte 
Carlo models. The first group comprises several versions of pythia: pythia 6 (tunes D6T, 
Zl.LEP LIZJ/ AMBTl, DW-Pro, and Pro-PTO provide very similar results), PYTHIA 8 (tunes 2, 
3, 4, and 5 are equivalent) and the recent PYTHIA 8 MBR tune llTSl . The second group includes 
MC generators based on the same Regge-Gribov phenomenology, but with different imple- 
mentations of model ingredients |19|: Phojet, as well as three MC programs commonly used 
in cosmic-rays physics, such as QGSJET 01 [20 1, QGSJET II (versions 03 and 04) ||2T1, SiBYLL 
(version 2.1) [22j and EPOS (version 1.99) Il23l . 

The Phojet and Sibyll models overestimate the observed cross sections by more than 20%, 
while the EPOS, QGSJET 11-03, PYTHIA 6, and PYTHIA 8 tunes provide predictions that are 
about 10% larger than the measured inelastic cross sections. QGSJET 01 and QGSJET 11-04 agree 
within one standard deviation with the data points. The pythia 8-MBR tune reproduces rather 
well the vertex-based measurements, while it overestimates the calorimeter-based result. 

A comparison of the trends in the data with the MC models is shown in Fig. |6| where the 
cross sections are now normalized to the (Tjnei value measured for events with >3 tracks. The 
dependence of CTj^ei on the nature of the final states relative to the results for >3 tracks, is well 
reproduced by most MC simulations. 

The TOTEM collaboration [2] has recently measured a total pp inelastic cross section of 

c^inel = 73.5+^ 9 mb. Although several Monte Carlo models such as EPOS, QGSJET 01, QGSJET 

II-4, PYTHIA 6, and pythia 8 reproduce this value (Fig.|5]l, only QGSJET 01 and QGSJET 11-04, 
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Figure 4: (a) Values of the inelastic pp cross section crjnei and (b) their associated goodness of fit 
;:\;^/NDOF, obtained for each of the fits in Fig-jsj as a function of the number of pileup vertices, 
in interactions with > 1 track with pj > 200 MeV/c and < 2.4. The line in (a) is the result of 
a fit to the 9 individual values of Oinei/ while the dashed line in (b) indicates ;\;^/NDOF = 1. 

and PYTHIA 8-MBR (but less so) are able to simultaneously reproduce the less inclusive CMS 
measurements. This observation suggests that most of the Monte Carlo models overestimate 
the contribution from high-mass diffraction to the total inelastic cross section, and underesti- 
mate the component at low mass. 



6 Summary 

The inelastic cross section in pp collisions at ^/s = 7 TeV has been measured using two methods 
that incorporate information either from central or from forward detectors of CMS. The results 
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Table 4: List of systematic sources and their effects on the value of the inelastic cross section 
measured using the vertex-counting method. The % changes are shown for the results of the 
^mel(>l track) measurement. The luminosity contributes an additional uncertainty of 4% to 
this measurement. 



Systematic source 


Uncertainty on cTj^ei 


Change in crjj^e\{>l track) 


Vertex reconstruction efficiency 


±1.4 mb 


±2.4% 


Longitudinal position of vertex 


±0.1 mb 


±0.2% 


Vertex quality 


±0.7 mb 


±1.3% 


Minimiun distance between vertices 


±0.1 mb 


±0.2% 


Transverse position of vertex 


±0.3 mb 


±0.6% 


Different sets of data 


±0.9 mb 


±1.6% 


Range of luminosity used in fit 


±0.2 mb 


±0.4% 


Reweighting MC track distribution 


±0.2 mb 


±0.4% 


Total (in quadrature) 


±2.0 mb 


±3.3% 



Table 5: a^j^gi values for interactions with >1, >2 and >3 charged particles, with their uncer- 
tainties from systematic sources of the method and from luminosity. 



Measurement 


Result 


c^mei(>l track) 
^^mei(>2 tracks) 
t^mei(>3 tracks) 


[58.7 ± 2.0 (syst.) ± 2.4 (lum.)] mb 
[57.2 ± 2.0 (syst.) ± 2.4 (lum.)] mb 
[55.4 ± 2.0 (syst.) ± 2.4 (lum.)] mb 



for the different choices of final states considered are: 

>^mei{^ > 5 X 10"*^) = [60.2 ± 0.2(stat.) ± 1.1 (syst.) 

±2.4 (lum.)]mb, 
(^inel{> 1 track) = [58.7 ± 2.0 (syst.) ± 2.4 (lum.)] mb, 
£rinei(> 2 tracks) = [57.2 ± 2.0 (syst.) ± 2.4 (lum.)] mb, 
£rinei(> 3 tracks) = [55.4 ± 2.0 (syst.) ± 2.4 (hxnx.)] mb, 

where each track must have pj > 200MeV/c and \rj\ < 2.4. The comparison of these results 
with the cross section expected from Monte Carlo models used in collider and cosmic-rays 
studies shows that Phojet and SiBYLL largely overestimate cTi^ei. The EPOS, QGSJET 11-03, 
PYTHIA 6, and PYTHIA 8 (except the MBR tune) programs predict values about 10% above 
the data, while Qgsjet 01, Qgsjet 11-04 agree well with the measurements. PYTHIA 8+MBR 
agrees well with the track-based measurements, but overestimates the prediction for (Tj^ei for 
^>5xl0~^. All models agree broadly with the dependence of the cross section on the criteria 
used to define the final states. 
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